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A detailed discussion is given of the disorder parameter for dual superconductivity of lattice gauge theories, introduced in a 
previous paper, and of its relation to other approaches. New lattice data are reported. Among other results, we find that the 
dual superconductivity of compact (7(1) is type II. 
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I. INTRODUCTION 

Dual superconductivity of the vacuum is an important 
phenomenon in gauge theories. 

- It produces confinement of electric charges via dual 
Meissner effect in the abelian case. 

- It is most likely the mechanism of colour confine- 
ment in QCD ||,| 

- It plays a key role in the structure of supersymmet- 
ric gauge theories [1). 

The simplest case is the compact U{1) gauge theory. 

With Wilson action this theory shows a phase transi- 
tion at (3c — 1.01, probably weak first order, from a phase 
at low (3 where electric charge is confined, to a phase of 
free photons [||. 

Confinement is detected by measuring the string ten- 
sion from the vacuum expectation value {vev ) of Wilson 
loops. The penetration depth of the electric field is fi- 
nite for P < Pc, indicating dual Meissner effect, and goes 
large at the deconfining transition Pc [Q. 

Monopoles are detected by their Dirac strings as units 
of 27r magnetic fiux through the plaquettes. Their num- 
ber density is not a disorder parameter for dual super- 
conductivity, in the same way as the number of electric 
charges is not for ordinary superconductivity. However, 
empirically, the number density of monopoles is larger in 
the confined phase, and drops to zero above pc 

A legitimate disorder parameter should vanish for sym- 
metry reasons in the deconfined phase, and be different 
from zero in the confining phase. Since dual supercon- 
ductivity is nothing but the spontaneous breaking of the 
U{1) symmetry related to the magnetic charge conserva 
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tion, the vev of any operator carrying non zero magnetic 
charge can be a disorder parameter. A non zero vev of 
such an operator would indeed indicate that vacuum has 
not a definite magnetic charge, i.e. that monopoles con- 
dense in it in the same way as Cooper pairs do in the 
ground state of ordinary superconductors. The concept 
of disorder parameter is known since long time in the 
comunity of field theory and statistical mechanics [^1). 
In the pionering numerical simulations in lattice gauge 
theories, however, emphasis was given to the density of 
monopoles as indicators of dual superconductivity . 

A rigorous proof was given in ref. |^ that monopoles 
condense at low /3's in lattice U{1) theory with Villain 
action ||]. The proof makes use of the specific form of 
the action and so did numerical attempts to extract a 
disorder parameter pc| ]. However, probably because of 
the mathematical language of the forms, which is not 
so familiar to physicists, nobody tried for long time to 
export the construction to the generic form of the action, 
or to non abelian gauge theories. Indeed, after abelian 
projection [ pT| , monopole condensation in non abelian 
gauge theories like QCD, always reduces to an effective 
U{1) with Dirac monopoles |ll|,|l^. Of course the U{1) 
effective action is unknown, and therefore a construction 
of the disorder parameter is needed, which can work with 
any variant action. 

Such a construction was given in ref and immedi- 
ately afterwards was used to demonstrate dual supercon- 
ductivity of non abelian theories [ p^ . 

This result prompted the exportation of the construc- 
tion of ref. |^,|l^ from Villain to generic action. 

In this paper we want to discuss in detail and improve 
the construction of ref |l^ (sect. 2), compare it to that 
of ref. 1^, (sect. 3), showing that they are equivalent, and 
present a number of numerical results for lattice U{1) 
with Wilson's action, (sect. 4). We will then compare our 
|]l3| way of detecting superconductivity by the quantity 
P = lii(/i), {fj,) being the disorder parameter, to direct 
determination of (/x) or of its effective potential (sect. 4). 

Besides confirming dual superconductivity of U{1) 
gauge theory in the confined phase we show that it is 
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II kind. 

The present discussion is also a useful basis to the 
treatment of the analogous problem in non abelian theo- 
ries, which will be presented elsewhere. 



II. THE DISORDER PARAMETER. 

The construction of ref of the creation operator of 
a monopole or antimonopole, is inspired by ref |^,^ and 
is based on the following simple idea. 

In the Schrodinger representation where the field A{x) 
is diagonal, a monopole of charge 27r- sitting in y is 
created by adding the corresponding vector potential 
^b{x~y) to A{x). 

This is nothing but a translation of A{x), which is 
generated by the conjugate momentum Tf{x) = E{x), the 
electric field operator. In the same way as 



we have 



with 



(1) 



(2) 



(3) 



(4) 



the commutator can be evaluated by use of the 

canonical commutation relations 



e'P^lx) = \x + a) 



\A{x,t) + -b{x-y)} = ^i\A{x,t)) 
e 



1 



1- / d'^xE{x, t)b{x- y) 



Ky,t) = exp 
The magnetic charge operator being 

d^xV f V A A{x, t) 



[E,{x,t),Aj{z,t)] = -i5^j8^{x- z) 



(5) 



givmg 



[Q(t), ^i{y, t)]^ - j d^xV ( V A 6(f - y) ) • ^i{y, t) 



= ^tx{y,t) / d?xV 



2e 

27r-M(2/,i) 



(6) 



In deriving Eq.(||) the Dirac string has been removed. 
A choice for b{x — y) can be 

o{x-y) = -- 



2 r(r — r ■ n^) 



(7) 



Alternative choices differ by a gauge transformation, b — s- 
h + which leaves the operator invariant if the Gauss 
law V • E' = is satisfied. 



On the lattice the building block of the theory is the 
link U^{n), which is an element of the gauge group. For 
U{1) Uf_,{ji) = e'^"'"^ and the plaquette, II^j^, which is 
the parallel transport along the elementary square in the 
plane fiv at the site n, is 



with 



9nu — ^ufiv — ^uSn — a eFni, 

a— >0 



The lattice version of the electric field is then 

e 

and a definition of the operator /i on the lattice can 
be 



(8) 
(9) 

(10) 



fJ-iy, no) = exp 



exp 



-^^5'(n-^)Imn°'(n, no) 



-/3 2^b''{n- y) sin(6loi(n, no)) 



(11) 



(3 = Here U{n) is the discretized version of the 

monopole field Eq. . The factor [3 in front of the ex- 
ponent comes from the factor ^ in the monopole charge 
times the normalization factor in Eq. ([lO|) . Usual Wick 
rotation to Euclidean region has been performed. 

The form (|ll|) was successfully used in ref. jl^ . 

A better definition of /i can be given, which coincides 
with Eq.(pT[) in practice, but automatically respects the 
compactness of the theory, in that it shifts the exponent 
of the links, and not the links themselves. In formulae: 

fJ-iy^mo) = exp|/3^ ( cos [61°' ™o) + b\y~ n)] - 

n 

- C08[e"\n,mo)])} (12) 



For small the definition ([l^ ) coincides with ([III). 

More generally if J2fj,^n S{0^y{n)) is the action, ^ will 
be defined as 

M(y,™o) = exp|/3^ [S'(6'*(^i,™o) + h{n~y)) - 

n 

- Sie°\n,mo))]] (13) 



and will tend to the expression (|llj) as the lattice spacing 
a go to zero, when the action tends to the continuum 
action. 

The prescription of excluding Dirac string on a lattice 
being either to locate the monopole at y between two 
neighbouring sites, or to eliminate in the sum the arrow 
of sites where b is singular, it is easy to verify that the 
definitions ^ and (^ give the same results from the 
practical point of view. 
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If the action is the Wilson's action ||15| 

J2 /^(cos(V)-l) (14) 

n,(/Ai>') 

then the vacuum expectation value of fi is given by 



or, making use of 
{fi{y, mo)) ^ 



n d^M(") 



n 



exp(S') /i 



(15) 



exp(5 + S') (16) 



where 5' is the exponent of Eq.(p^. 

Adding 5' simply amounts to modify the (0, i) plaquet- 
tes on the time slice np, by addition of bi to 0oi 



(17) 



S + S' = J2 J2 /3(cos(%(n) -1) + 
+ J2 P (cosiOo^in) ~ 1) + 
+ ^ (cos(6'oi(^, mo) + b\rn - n)) - l) 



If a number of monopoles and antimonopoles are cre- 
ated at time no, bi should be the sum of the correspond- 
ing vector potentials. The generic correlation function 
{fj,{xi) . . . fi{xn)) is defined as (fi) in Eq.(p^, with the 
change from S to S + S' extended to all the time slices 
where monopoles or antimonopoles are created. 

So for example the correlation function where a 
monopole is created in y = at t = and destroyed 
at i = no is given by 

{^iiy, 0) li{y, mo)) = ^ / exp(5 + S'^^) 



S + S' differs from S by the replacement 

6lo,(n, 0) -> %(n, 0)+b,{n~y) at i = 
9Q,{n,mo) ^ 9oi{n,mo) - bi{n - y) att = mo (18) 

Monopole condensation can be detected from the asymp- 
totic value of (/i(y, 0) fl{y, no)). Indeed as no grows large, 
by cluster property 



(M(y,0)Ai(27, »7io)) ^ C exp(-moM) + {^lY 



(19) 



Notice that (/i) — (p) by C invariance, and the position 
y is irrelevant by translation invariance. M is the mass 
of the lowest state with monopole charge q in units of 
inverse lattice spacing. 

To visualize that fi really creates a monopole at t — 
consider again the change it produces according to 
Eq.(|l8|). Since 



00,: (n, 0) = e,{n, 1) - 0oin, 0) - 0o{n + I 0) + 0o(n, 0) 

(20) 

the change ( |l8|) of 9oi can be considered as a shift 

e,{ri,l)^e,{ri,l)-b,{ri-y) (21) 
A change of variables 



9i{n,l) - bi{n - y) 



(22) 



in the Feynman integral (pl|), which leaves the mea- 
sure invariant, brings back the plaquette O^i to its un 
perturbed form. However the change of variables (|2~ 
changes the (i, j) plaquette at no = 1 as follows 



MA) 



,(n, 1) + A»&j(n -~y)~ Aj-6,(n - y) (23) 



This means that at no = 1 the magnetic field of a 
monopole located aX n — y is added to the original con- 
figuration. The change of variables ( pl| ) also affects the 
plaquette 9Qi{n, 2), and amounts to the shift 



^o^{n,2)^9o^{n,2)~b,{n~y) 



(24) 



Again a change of variables Oi{n, 2) Oi{n, 2) — bi{n — y) 
restores 9oi{n, 2) to the initial form at the price of adding 
a monopole at time t = 2, and of producing a shift in the 
form ( p^ ) on 6oi{n,3). This procedure can be iterated. 
At t = niQ this procedure ends, because bi cancels with 
the shift of opposite sign corresponding to the creation 
of the antimonopole. 

Thus the correlator {fi{y,0)fl{y,'mo)) simply consists 
in having a monopole propagating in time, from to no. 

The construction above simply generalizes to more 
complicated forms of the action, where Wilson loops 
other than plaquettes enter. 



III. COMPARISON WITH OTHER 
APPROACHES 

In this section we want to discuss the relation of our 
approach to that of ref. 

In the language of ref. 1 1| 0^ (n) is a 1 form associated 
to the links and d6' is the two form associated to the 
plaquettes, or the field strength tensor. 

In this language the partition function is 



Jv[9] $^(d0) 



For Wilson's action 



$^ = exp(/3^ (cos(d6i) - 1) 

plaq 



(25) 



(26) 



For Villain's action 
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plaq 



(27) 



To define a disorder operator (fi) the action is modified 
by adding a two form X to d6. We define 



and 



Z{X) = J V[e] $/3(d6i + X) 
Z{X) 



m 



(28) 



(29) 



Any change of X of the form X ^ X -\- dA leaves (/i) 
invariant, in that dA corresponds to a shift of 6* to 6* + A 
which is reabsorbed by a change of the (periodic) inte- 
gration variables. 

Since a generic X can be written as (Hodge decompo- 
sition): 



X = da + S^dX 
A 



(30) 



the above invariance implies that (/i) only depends on 
dX. 

dX is a 3- form [dXj^^Q and its dual *dA is a 1 form, 
which is a magnetic cm'rent, since A is a field strength. 
Explicitely 



dApm — — {daXf^u + dfj,Xua + dyXaf^) (31) 



and 



J, 



6 



(32) 



The magnetic current ( p2D is identically conserved. In 
the language of forms 



6 J'' = 



(33) 



The magnetic charge density which describes the creation 
of a monopole of charge 2Trq in the site y at time y^, and 
its destruction at time y"^ is 

Jo*^(f, x°) = 27rqS^x~y) (0(.t° - y°) ~ 0{x' ~ y'")) 

(34) 

Since the current is conserved 

VJ^^ = -Ao J(?^ = 

= -27rqdHx-y) {Six^ - 2/°) - Six' - y'°)) (35) 

A solution of Eq. ( ^ is 

1 x-y 



J^'{x,X°)^2^q 



An 



x~y\ 



[<5(x«-y°)-5(a;°-y'0)] 
(36) 



The corresponding X is then 

X^5^J^" 
A 



(37) 



The correlation function of a monopole antimonopole will 
then be: 



Z{0) 



(38) 



This is the construction of ref. [g| . 

Notice that Z{X) is periodic in X (with period 2tt) 
since the action is compact. In fact Z only depend on dA, 
and is periodic also in dX with the same period. This 
can be rigorously proved by going to Fourier transform: 

Z{dX + 2TTn) = Z{dX) (39) 

Consider now a one form on the dual lattice, with 
support on a line. If Sfl — the support must be a closed 
line. If il is integer valued in units of 2t: the change 

dA = * j^^ dA = + n 



leaves Z invariant. 

In the notation of ref. 
Jq^ by ~2TTqL0 and 

dA = 2Tiq{B 



J is denoted by 2TiqB and 



(40) 



Any X with the same dA, will give the same correla- 
tion function (|3^). The construction presented in sect. 2 
corresponds to the choice 



Xr 



h{x) [5{x'-y')-5{x'^-y'')\ 



X,, = 
or, in the dual language 
(*a')o, = 

(*a'),, = e,,khk{x) [<5(xO - yO) - J(a;0 - y'^)] 
and 



*dA, 



5UX 



p 



Ap(*A), 



(41) 

(42) 
(43) 



Explicitely 

(5(*a')o = 



5(*a'), = -^Afe(*A) 



ki 



and by Eq.'s (|T|) and (0) 

5(*a'), = 27rg— ^^^^ 



A I- ^3 ('5(2;o - yo) - 

47r |a; — y] 
2nq6{xi - yi)5{x2 - y2)0{x3 
{6{xo - 2/o) - S{xo - y'„)) 



6{xo-y'o)) (44) 

- ys) ■ 



Our *dA differs from *dA (|o|) of ref. [| by a 1 form 
integer valued in units of 27r, with support on a closed 
line. Therefore our correlator coincides with that of ref. 

not only for Villain action, but for generic form of 
the action. 

This section is a cultivated way of presenting the ar- 
gument already given at the end of last section. 
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IV. NUMERICAL RESULTS FOR THE 
DISORDER PARAMETER. 

As discussed in sect. 2, we measure the correlation 
function 

V{x") = {fi{x, x"^), p{x, 0)) ~ Ae-^'^° + {fif (45) 

The aim is to extract (/i)^, which wiU signal dual su- 
perconductivity, and M which is the lowest mass in the 
sector of magnetically charged excitations. 

A direct determination of V can be done, as we will 
discuss below, but is rather noisy from numerical point 
of view. The reason for this is that V 



P = i Jve cxp{S + S') 



(46) 



is the average of exp(5"), S' being the modification of 
the action on the time slices t = and t = x'^, and 
S' fluctuates roughly like the square root of the spatial 
volume. 

A way to go around this difficulty is to measure, in- 
stead of V the quantity |12 



p{x, X, 0) = — In I? 



dp 



(47) 



At large distance (a;° — > oo) 

p^^2^\n{^i) (48) 

and since p{f3 0) = 1 (p.) can be reconstructed as 

1 



(a*> = exp 



From Eq.(p 



{S)s-{S + S')s+s' 



(49) 



(50) 



The definition of p is analogous to the definition of the 
internal energy in terms of the partition function in sta- 
tistical mechanics, p is now a well defined quantity and 
easy to measure, and, as we shall see, can give all the 
information needed to detect dual superconductivity. 

We have made simulations on a 6'^ x 12, 8'^ x 16 and 
lO'^ X 20 lattices putting the time axis along the long 
edge of the lattice. A typical behaviour of p versus x° 
is shown in Fig.l, for a 8"^ x 16 lattice, showing that an 
asymptotic value is reached by p as a function of . The 
mass M of the exponential in Eq.(p5|) can be estimated 
and is typically (2 — 3) /a 
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Fig.l Monopole antimonopole correlation in time. 
(Lattice 8^ x 16) 

We will come back again to this point in the following. 
The quantity poo as a function of f3 is plotted in Fig. 2. 
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1.2 



Fig. 2 Poo as a function of f3. The negative peak signals 
the phase transition. (Lattice 8'^ x 16) 

For all of our lattices sizes Poo is negative and sha rply 
decreases approaching Pc- This corresponds, by Eq.(|49|) 
to a behaviour of (p) which slowly decreases from the 
value (p) — 1 at (3 — 0, and has a sharp drop at Pc- 

To better analyse this behaviour we compare it for the 
three lattice sizes under study. For P < Pc below the 
negative peak, p increases with L, showing that as L — > 
cxD, (/i) reaches a finite, nonzero value. Magnetic U{1) 
is therefore spontaneously broken, and for P < Pc the 
system is a dual superconductor (fig. 3). 

For /3 ~ /3c we know that the typical correlation length 
of the system goes large. There is evidence that the tran- 
sition is weak first order ||l6|| , with some controversy [pT| . 

The correlation length ^ goes large as P approaches 
Pc in a range of /3's and eventually stops growing before 
reaching it. 
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Fig.3 Poo versus 1/i for /31.009. 

This means that, in the neighbourhood of (Sc 



M = M 



(51) 



If the transition were second order a critical index v 
would exist such that 



(52) 



In our case some effective index v could anyhow exist, 
describing a behaviour of ^ of the form ( |5^ ) in the above 
mentioned range of /3's. Then ^/L can be traded with 
L^/''{l3c - /3) and a finite size scaling behaviour results 



implying for p = jjs I^(A') ^ scaling behaviour 



-l/u 



(53) 



(54) 



Eq.(p4|) allows a determination of v and /3c, togheter with 
a determination of the exponent 6 by which (/i) tends to 
zero at /3c in the infinite volume limit. 

The quality of the scaling is shown in Fig. 4. points 
corresponding to different lattice sizes follow the same 
universal curve only for the appropriate values of /3c and 

Eq.(|55|). If f3c or i/ are changed by one standard de- 
viation from the values of Eq.(p5|) points from different 
lattices start splitting apart from each other. 




10.0 20.0 

(Pc-P) l"' 

Fig. 4 Finite size scaling. pL^^^ is plotted 
vs. (/3c-/?)lV-. 



30.0 



A best square fit gives 



Pc = 1.01160(5) 
u = 0.29(2) 



(55) 



The value (55) of fic is consistent with determination 
based on completely different methods jl6j. If ^ ^ 

(/3c - /3)' 



(56) 



L^/- L1/-(/3c-/3) 
An estimate for 5 from the behaviour in fig. 4 is: 
5= 1.1 ±0.2 



(57) 



In the region /3 — > cxd p can be computed in the weak 
coupling approximation p3|. The result is 



-5.05 -i + 4.771 



(58) 



giving p — > — cxD or (/i) = in the infinite volume limit, in 



agreement with general arguments | 
the disorder parameter vanishes in t 
if boundary conditions are not free. 



: only as 1^ oo 
le disordered phase. 



The mass of the monopole in Eq. ( 45 ) should scale prop- 
erly in the limit (3 ^ (3c but we have large errors and this 
behaviour is not clearly visible (fig. 5). 

In order to determine if the superconductor is first kind 
or second kind we have also measured the penetration 
depth l/ruA of the electric field on the lines of ref. 
A constant electric field parallel to the space boundary 
of the lattice is put on a face of the space lattice and its 
value is determined inside the bulk as a function of the 
distance from the boundary. An exponential behaviour 
is found, with a penetration depth which properly scales 
by approaching the critical point, consistently with the 
effective critical index. 

The corresponding mass is shown in fig. 5 together with 
the mass extracted from the correlation length, Eq. (^ . 
It appears clearly that M > 2mA, indicating that the 
superconductor is second kind. This same problem has 
been approached by looking at the Abrikosov flux tubes 
generated by propagating charges. The idea is to com- 
pare the dependence of the electric field inside the tube on 
the transverse distance x± from the center of the tube, 
with what is expected from London equations. Their 
result is that the system seems to be is at the border 
between first and second kind jl^ . The method is ingen- 
uos. However derivatives are approximated by finite dif- 
ferences, the penetration depth being a few lattice spac- 
ings (2-3), and this can produce systematic errors. Our 
method would give a more precise determination if we 
were able to determine better the mass AI of Eq. (^5|) . 
The question deserves further study. 



^In ref. Q the field was called magnetic. 
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Fig. 5 Mass of the monopole M (squares), and mass of the 
dual photon to (circles) vs. /?. 

Finally we want to comment on the possibility of de- 
termining numerically (/i) directly and not through the 
measurement of p. As we have seen this is not strictly 
necessary, since p gives complete information about the 
phase transition. However the problem has some interest 
by itself. 

The definition of (/.t) is (e'^"^ ), the average being per- 
formed with the weig ht VOef^^ /Z. S' is itself a random 
variable in this ensemble which has some average value 
{S') with a width a = .JJS^^y^n^. 

A general theorem of probability theory states that if 
a random variable is distributed with a probability law 
p{x), with J p{x)dx — 1, then its average Xn = ^ Tlk^ 
is distributed as a gaussian for large n if and only if |19| ] 



lim 



p{x)dx 



x\>X 



x^p{x)dx 



(59) 



l\x\<X 

If Eq.dl) holds, then 



{Xn) -> (x) = / Xp{x)dx 



(60) 



and the width of the distribution is in this limit 

a 

with 



theorem of central limit. In fact a much slower decrease 
would be enough. 

If, for the sake of the argument, we assume that 
is gaussian, then we easily compute, by use of Eq. ( |6l|) 



(a*) = Mexp(-|-) 
(T^ = /2exp((T^) 



(62) 



Eq.'s(62) show why a direct determination of (/i) is af- 
fected by wild fluctuations: the width is indeed bigger 
than the value of (^) itself. The exponential dependence 
on S' strongly distorts the distribution when going from 
S' to /i. 

The histogram of the values of p, is related to the con- 
strained potential by the relation [fo 20 



exp (-y(«')) = y \DB\ exp{f]S)d{p - $) 

F($) has a minimum at (/3S") -I- 

If instead we construct the histogram of /3S" itself, the 

minimum will appear at {PS') which is displaced by -j- 
with respect to the real minimum. 

The problem is that the histogram in p is exponen- 
tially large to fill adequately, since p fluctuates on an 
exponential scale (typical values of p on a configuration 
for a reasonable lattice size range from 10^^" to 0). A 
histogram of log p, i.e. of f3S' is easier to compute. 

However to go back to the distribution in p, i.e. to 
compute (p) and cr^, we must know the distribution 7r(y) 
with great precision. In the gaussian approximation the 
solution is given by Eq.(|62|). A cluster expansion can 
be attempted, to evaluate non gaussian effects, but the 
problem is only shifted. Higher cumulants of 7r(?/) are 
more and more noisy to determine numerically, and the 
computer time needed becomes comparable to the one 
needed for the direct determination of (p). 

Finally a finite size scaling analysis would be needed, 
analogous to what we did in sect. 4. 

This is to justify why we used p to extract information 
on the phase transition, instead of (p) itself, or of its 
effective potential. 

The problem is currently under further study. 
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{x~ 



{X 



If we denote by y the variable (3S' — {(3S') and by 7r(y) 
its probability distribution, then the variable p 



p — e'xjp{pS') — pexp(y) 
will be distributed as 



(a2 = expiiPS'))) 



p{p) = TT ( In ( - 



dln( ^ 



(61) 



If TT decreases as exp(— y^/2cry) as y ^ oo then the prob- 
ability distribution (plj) obeys the hypotesis (poh of the 
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